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1 Introduction 

The study of hard exclusive reactions in the Bjorken limit is crucial to obtain 
information on parton dynamical correlations in the nucleon. The recent data 
from DESY ep collider HERA on exclusive diffractive virtual Compton process 
|1I2I3 | (DVCS) at large becomes an important source to study the partons, 
in particular gluon, inside the proton for non-forward kinematics and its relation 
with the forward one. A considerable interest of the DVCS process comes from the 
particular access it gives to these generalized parton distributions (GDP) through 
the interference term with the Bethe-Heitler process. On the other hand, recently 
the color dipole formalism has been provided a simultaneous description of photon 
induced process. The inclusive deep inelastic reaction and the photon diffractive 
dissociation has been successfully described and the study of other exclusive process 
as DVCS is an important test of the color dipole approach. The work reported here, 
summarizing the studies in Refs. j4l5j . applies the successful saturation model [H] 
to the DVCS process. The model interpolates between the small and large dipole 
configurations and has its parameters obtained from an adjust to small x HERA 
data. Moreover, its QCD evolution has been recently computed '7^, which improves 
the high data description. 



2 DVCS cross section in the color dipole picture 

Based on the color dipole framework, the DVCS process can be seen as a succession 
in time of three factorisable subprocesses: i) the photon fluctuates in a quark- 
antiquark pair, ii) this color dipole interacts with the proton target, iii) the quark 
pair annihilates in a real photon. As usual, the kinematic variables are the c.m.s. 
energy squared s = = (p + 9)^1 where p and q are the proton and the photon 
momenta respectively, the photon virtuality squared = —q^ and the Bjorken 
scale X — Q^/{W^p + Q^). The DVCS imaginary part of the amplitude at zero 
momentum transfer reads as, 

ImA{s,Q\t^O)^ dz J d^r^*r{z, r, QI=Q^)^t{z, r, Ql ^ 0) <7d^p{^, r^) , 
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Figure 1. (a) Comparison between the profile function W{r,Q^) as a function of the dipole size 
r for DVCS and inclusive DIS processes at = 10 GeV^ . (b) The ratio of the imaginary parts 
of the DIS and DVCS amplitudes as a function of Q'^ at W = 75 GeV. 



where adip{x, r^)^ is the dipole cross section, which depends on the scahng vari- 
able X = ^^a^^g^) and dipole size, r. The product of the photon wavefunctions 
(transverse polarization) is given by, 

^ / 

where £12= ^{^ ^ z)Q\ 2 + rrij and f{z) — [z^ + (1 — z)'^]. The quark mass, m/, 
plays the role of a regulator as — *■ 0. The relative contributions from dipoles of 
different sizes can be analyzed with the weight (profile) function, 

W{r, Q") = / dzr ^'?,(z, r, Qj = Q^) ^t(2, r, QI = 0) a^p (i, r^) . 
Jo 

In the DIS process the contribution of large dipole configurations is observed 
to diminish in a sizable way as the virtualities increase. However, as shown in Fig. 
(^a), the DVCS profile selects larger dipole sizes in contrast to the inclusive DIS 
case, even at relatively large (similar result has been obtained in Rcf. 8 ). The 
calculation was performed using the saturation model "H* for both three and four- 
flavor analysis. The inclusion of the charm content gives a lower normalization for 
the profile and by consequence for the total cross section. The impact of the charm 
is smaller in the inclusive DIS case than in DVCS process, confirming that DVCS 
is more sensitive to the non-perturbative (soft) content of the scattering process. 
Moreover, in order to estimate the importance of the skewing effect, we calculate 
the ratio between the imaginary parts of the forward t — Q amplitudes for DIS and 
DVCS, R — 2m DIS/ Im DVCS. As shown in Fig. iQ-b) our result presents values 
slightly above those from an aligned jet model analysis in Ref. [5] and below those 
from the dipole analysis in Rcf. ~E . 

The final expression for the DVCS cross section is written as, 

aij p^^p) = ^ (l + p), 
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Figure 2. The photon level DVCS cross section as a function of virtuality at W = 75 GeV. 
Data are compared to the present prediction, with and without (BGBK) QCD evolution (see text). 



where B is the t slope parameter and comes from a simple exponential parameter- 
ization and p, the ratio between the real to imaginary part, is computed according 
to the dispersion relations (more details in Ref. 0]). 

3 Cross section comparison to Data 

In what follows we present the results for the phenomenological saturation models 
and also its version considering QCD evolution [7] (labeled BGBK), giving the 
dipole cross section gluon dependent. In Fig. [21(from the HI conference paper [5]) 
one compares their prediction with the existing measurement from HI and ZEUS 
Collaborations |1I2I3) . The color dipole prediction from Donnachie-Dosh is 
also presented for the sake of comparison. As the B value has never been measured 
for DVCS, the normalization of the theoretical prediction is basically free (usually 
values of 5 < B < 9 GeV~^ are considered, and the fixed value of i? = 7 GeV~^ is 
chosen on the figure). 

Although there is a little difference of normalization between HI and ZEUS 
measurements, which makes difficult to set an overall B value for all measurements, 
the behavior on Q"^ and W (see Ref. [3]) is well reproduced. On the other hand, 
for Q^>AO GeV^ our prediction still underestimates the experimental data. This 
change of behavior in the shape can indicate two situations: (a) the B slope 
would diminish as increasing virtualities or; (b) some additional effect should appear 
at higher Q^. In order to investigate the first hypothesis, we compute cross section 
using a dependent slope: B{Q^) = 8[1 - 0.151n(QV2)] GeV^ (see Ref. 
for details). Concerning the second hypothesis, we have investigated two options: 
QCD evolution (using BGBK model) and skwedness effects. For the skewedness 
corrections, the ratio of off-forward to forward gluon distribution are are given 

iT^i 2A-i-3 r ) 

explicitly by JI], Rg (Q^) = r (a+4) ' where A is the effective power on energy 

of the scattering amplitude. For our purpose the amplitude is multiplied by Rg, in 
order to estimate the size of the skewedness effects. 

To compare the dependence, we normalize all models to describe the ZEUS 
data point at the lowest value. Further, wc plot the ratio of each model to 
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Figure 3. The ratio MODELS/SAT-MOD as a function of (see text for details). 



our baseline model SAT-MOD as a function of . Such a procedure allows a 
dependence comparison independently of the normalization effect. These ratios are 
shown in Fig. O where the points (triangles-up) are the ratio of the ZEUS data to 
SAT-MOD including the error bars for the statistical (inner) and sum in quadrature 
of statistical and systematic (outer) uncertainties. 

We verify that several models can account for the measured dependence, 
which are not distinguishable with the present experimental precision. If the change 
in normalization is small for the inclusion of a dependence in B, the effect 
is of the order of 12% for BGBK with respect to the basic SAT-MOD and of 
40% for the skewedness effect (SKEW) and still larger when the different effect 
are combined (60% for BGBK-I-SKEW). Therefore, these issues show clearly the 
importance of a measurement of the slope B. Such a measurement would already 
allow to discriminate among the different theoretical predictions with an amount 
of data comparable to the present ZEUS measurement. 
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